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Abstract 

t-H _ 

We report on our reexamination of the radiative corrections to 7^2 and K12 decays. 
We perform a matching calculation, using a specific model with vector meson dominance 
in the long distance part and the parton model in the short distance part. By considering 
^> I the dependence on the matching scale and on the hadronic parameters, and by comparing 

£f} . with model independent estimates, we scrutinize the model dependence of the results. For 

the pseudoscalar meson decay constants, we extract the values / ff = (92.1 ± 0.3) MeV and 
fx — (112.4 ± 0.9) MeV. For the ratios R n and Rk of the electronic and muonic decay 
modes, we predict R^ = (1.2354 ± 0.0002) • 10~ 4 and R K = (2.472 ± 0.001) • 10" 5 . 
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< >i-j ! 1 — Introduction 

X: 

1.1 What is interesting about tt^, Ki 2 decays? 

Interest in 7T/2 and K\2 comes from two different sources. On the one hand, the decay widths 
of the pion and the kaon into muon and neutrino are used to extract the decay constants f n 
and f K , 

r(7T -> flUn) — ► f n 
T(K -»■ nu^) — ► fx 

which are important input parameters for chiral perturbation theory. And so the question 
arises how 0(a) radiative corrections affect these these parameters. 

On the other hand, the ratios and Rk of the electronic and muonic decay modes 

r(7T -» ev e ) T(K -» ev e ) 

can be predicted by theory. What makes R^ and Rk very interesting is that hadronic uncer- 
tainties cancel to a very large extent in these ratios and that the electronic decay channels 
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are strongly helicity suppressed in the standard model due to the V — A structure of the weak 
interaction. Thus the ratios R n and Rk are very sensitive to non standard model effects 
(such as multi-Higges, leptoquarks or other Q), which might induce an effective pseudoscalar 
hadronic weak current. And so 7172 and K12 decays allow for low-energy precision tests of the 
standard model, if we are able to understand the radiative corrections SRqed'- 

Rn = — X ItoI 1 + SRqED) ( 1 + 5Rnon-SM ) (2) 

The experimental precision tends to be better for 7172 decays. However, non standard model 
effects SRnvn-sM > might be enhanced in Rk by a factor of mx/m-ir, and therefore K[ 2 are of 
comparable interest in testing the standard model. 



1.2 Current Experimental Situation 

For R n , there are recent results (from 1992/93), obtained by TRIUMF and PSl(2j, resulting 
inpl a relative precision of 

^ = 3 x 10- 3 (3) 

The precsion of Rk as quoted by the 1994 particle data book 

AR K A irt _ 2 



Rk 



4 x W~ z (4) 



is one order of magnitude less, but one should note that this number is based on 1972-76 
data. It would be very interesting indeed to receive some new experimental information on 
this important quantity Rk- 



1.3 Theory 

There is a vast literature on the radiative corrections to 7172 decays. In 1958/59, Berman and 
KinoshitaQ calculated the radiative corrections to 7172 decays using a model with an effective 
pointlike pion. There have been numerous later attempts to improve this simple model in order 
to take into account hadronic structure effects. Important milestones in the understanding of 
these effects were the 1973 paper by Terent'ev^], who derived a general theorem on the leading 
hadronic structure dependent correction, and the 1976 paper by Marciano and Sirlin||, who 
proved that the leading logarithm of the radiative correction is not affected by the strong 
interaction, if one considers the inclusive decay rate T(7r — > fxu^) + r(7r — ► fiv^) integrated 
over all photons. In 1977, Goldman and Wilson^] considered various models for the strong 
interaction effects in the context of gauge theories of the weak interaction. In 1990, Holstein|8|] 
was the first to point to the importance of short distance radiative corrections in the extraction 
of f n , and in 1993, Marciano and Sirlin reconsidered 7172 decays from a modern point of view, 
collection all the leading model-independent corrections and giving rough order of magnitude 
estimates for the remaining model dependent corrections. 

Last year, we have performed an improved matching calculation]^] of 7172 and K\2 decays. 
In this calculation we include vector mesons as explicit degrees of freedom (instead of matching 
simply the point pion/kaon with the short distance part), calculate the full short distance 
correction (instead of including only the leading logarithm) and we scrutinize the size of 
the model-dependent corrections by explicitely calculating them within a reasonable model 
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(instead of just roughly estimating their order of magnitude). This calculation allows to 
obtain improved values both for the pseudoscalar decay constants f n and fx, and for the 
ratios R n and Rk, with a clear controll over the (small) model-dependent contributions. 

2 — Outline of the Calculation 

2.1 General Considerations 

To obtain the 0(a) radiative corrections to tt — > fJ,Un, we have to evaluate Feynman diagrams 
with one photon loop, where the photon can be contracted twice, once or not at all to the 
leptonic side. Because of the infra-red divergences, we have to add to this the tree diagrams 
for 7r — > iiVfTi, where the photon couples either to the leptonic or to the hadronic side. The 
essential task is then to understand the coupling of a (real or virtual) photon to the hadrons. 

In the treatment of the strong interaction, one can distinguish three different regimes. 
In the low energy regime, where momentum transfers are small compared with a typical 
hadronic scale of about 1 GeV, there are the low energy theorems of QCD which give model 
independent predictions. In the high energy regime, where the virtual photon momentum is 
large on a hadronic scale, we can use the parton model and calculate the QED corrections 
to the lv\ud operator. In the intermediate energy regime of momentum transfers of the order 
of 1 GeV, however, there is no way to obtain model independent predictions. This energy 
regime is dominated by the physics of vector meson resonances and can be described by 
phenomenological models only. 

And so in order to evaluate the loop diagrams, one separates the integration over the 
loop momentum into a long distance part h\ = • • • fi 2 ut and the short distance part k\ = 
f^cut ' ' ' m z> where [i cu t is a hadronic scale of about jjL cut ~ 1 GeV. 

Then one can either choose to neglect the problem of the vector meson energy regime and 
match the low energy regime and the short distances directly at a scale of about fj, cu t = rrvi, 
as has been done by Holstein|8| and by Marciano and Sirlin(T^]. A matching scale of m 2 , 
however, is somewhat too large for the low energy regime, on the one hand, becaue for 
P 2 — ► mi, all the neglected higher orders 0(P 6 ), 0(P 8 ), ...become large. On the other 
hand, m 2 is somewhat too small for the short distance part, where, after all, one is using 
asymptotically free quarks. 

Therefore we used the alternative approach of including the vector mesons as explicit 
degrees of freedom in the long distance part. This allows to push up the matching scale up 
to 

i^cut — 1 • • • 2 GeV, rendering the short distance part more reliably. 

However, this unavoidably introduces model-dependence. Thus we then scrutinize the size 
of the model dependence in different ways. Firstly we vary the matching scale from [i cu t = 
0.75- --3GeV and the parameters of the hadronics within their experimental uncertainties. 
Secondly, we determine the full size of the contribution to the loop integral from intermediate 
scales (0.5 • • • 3 GeV), which is the regime where the model dependence is large. Thirdly, we 
compare our result with the leading model independent contributions and then take the full 
size of the model dependent contribution as a measure of the uncertainty. 

2.2 Long Distance Part 

For the long distance part, we start from a model with an effective point meson model, which 
should be a good approximation for very small momentum transfers. These point meson 
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amplitudes are then modified along the lines of vector meson dominance in order to account 
explicitely for the vector meson degrees of freedom. Consider, for example, the coupling 
V(p,p')^ of the photon to a pion line. In the point meson approximation, this coupling is 
given by scalar QED as 

V(P,P%aint = i <P + P'y ( 5 ) 
where p and p' are the ingoing and outgoing pion momenta. But actually this coupling defines 
the pion electromagnetic form factor F n via 

V(p,p'T phy s = ieF n (k 2 )(p + pr (6) 

where V£ h denotes the physical photon-pion-pion coupling and k = p' — p is the photon 
momentum. F w (k 2 ) is known experimentally rather well in the relevant regime of k 2 up to 1 
or 2 GeV. Thus in the point pion amplitudes, we replace V(p,p > )^ oint by V(p,p')^ hys , using the 
parametrization of F n determined injn]]. This modification determines in turn the appropriate 
modification of the pion-photon-W boson seagull coupling, viz. by the requirement of gauge 
invariance. In the case of the kaon, the scalar QED point kaon-kaon-photon coupling has to 
be replaced by a coherent superposition of p, u and <3? vector meson dominance. In addition 
to these modified point meson diagrams, there are loops which correspond to the so-called 
structure dependent (SD) radiation in the radiative decays ir(K) — ► Ivff, where the emitted 
photon becomes contracted with the leptonic side. 



2.3 Short Distance Corrections 

For photons with large virtuality, \k 2 \ > (1 ... 2 GeV) 2 , we can use the parton model. Thus the 
first step is to calculate the one- loop photonic corrections 5 A to Aq = ^7= [u v ^^-Ui] [^7^7- u M ] . 

Neglecting all masses except for mi and Hcuti we obtain the leading logarithm of the 
radiative correction 

XA a 1 ( 2i m Z 2 , m z\ . 

77 rn l H-cut \ rn l H-cut / 

where the dots indicate corrections which are not leading in the limit m 2 z — ► 00. 

In the radiative corrections to the individual decay rates M — » li/ T (j), this leading loga- 
rithm clearly dominates, and so for the extraction of f n and fx its consideration is sufficient. 
It depends, however, only very little on the lepton mass and thus cancels almost completely in 
the ratios R n and Rk- In view of the very high precision of the theoretical predicton for these 
ratios, we therefore go beyond the leading logarithm and calculate the full 5 A. This full result 
for 5 A is firstly no longer proportional to the Born amplitude Aq and secondly it depends 
on the relative momentum of the two quarks. Therefore we firstly project onto the J p = 0~ 
component of the two quarks, and secondly we integrate over the relative momentum up of 
the quarks in the infinite momentum frame (u = — 1 . . . + 1). This leads to a result 

( SRm ) h +a- , = TT \ du $M{u)r M (u) (8) 

v / short dist £JM J — I 

Here $m(u) denotes the parton distribution function, whereas tm(u) is calculated from the 
projected short distance diagrams. We find that r n (u) and rx(u) depend only very little on 
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u, and so we can approximate them by their values at u = 0, where the wave function is 
presumably peaked: 



(SR M ) « r M (0) -f- f +1 du<S> M (u) = r M (0) (9) 

\ / short dist 6JM ■'— 1 



M J-l 
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where the last equation follows from a sum rule[p^|. 

In matter of fact, we find that R n and Rk are also dominated by the leading logarithm, 



(SR M ) ^A'^ do) 

V /short dist 7T - fi^ut ^cut 



2.4 Treatment of Real Photons 

Because of the infra-red divergences, one has to add some or all real photons to the virtual 
corrections and to consider inclusive decay rates 

T(M^lvi)+T(M^lva) 

Whether one includes only soft photons or all real radiation is a matter of convention. Our 
convention, which appears to be in accord with the one used by experimentalists, is to include 
all photons in the pion decays r(-7r — > £1^(7)). In the case of the kaon decays T(K — ► ^(7)), 
we include all the internal bremsstrahlung (IB) photons, but exclude completely the structure 
dependent (SD) radiation. 

For K — > fiv^jpf this SD radiation is completely negligible, but for the electronic mode it 
is extremely large, Tsd(K — > eUe'y) « To(K — ► eu e ), where Tq denotes the Born amplitude. 
Therefore it is usefull to consider the SD radiation as a separate decay mode and not to 
include it in the radiative correction to K — * ev e . 

Of course it is not possible to tell definitely whether a radiated photon is due to internal 
bremsstrahlung or to structure dependent radiation. However, if suitable experimental cuts 
are used, which put a small upper limit onto the photon energy, the measured rate of K — > 
ev e (-f) will include only a very small SD background, and only very little of the IB part will 
have been discarded. Using the predicted differential distributions JOg, [H| ], the SD background 
can be subtracted and the missing IB part added. Because of the smallness of this correction, 
it does not give rise to any important uncertainties. 



3 — Numerical Results 

3.1 The pseudoscalar decay constants 

Adding up the long and short distance corrections, we obtain the full radiative correction. 
The numerical result depends on the choice of the matching scale [i cu t and on the parameters 
of the hadronics, which are known only with limited experimental precision. 

In Fig. m we display the correction to the decay rate T(K — » ^^(7)) in variation with 
Mcut) using three different choices for the hadronic parameters. The solid line (I) corresponds 
to the central values of the hadronic parameters, whereas the dashed (II) and the dotted 
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Figure 1: Radiative correction to T{K — ► (iV/i), using different choices for the hadronic 
parameters: Standard choice (I, solid) and variations (II, dashed and III, dotted) 
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<0 I r — I 



^ CU i [GeV] 



(III) lines are obtained by varying these parameters within reasonable ranges. From this, we 
obtain 

— {K -» M^u(7)) = (1-23 ± 0.13 ± 0.02)% + 0(a 2 ) + 0(aa s ) (11) 
To 

The central value 1.23% of the radiative correction of order 0(a) has been obtained with 
fJ'cut = 1.5 GeV and the central values of the hadronics. The first error quoted, ±0.13% is the 
matching uncertainty, obtained from varying the matching scale up and down by a factor of 
two, fj, cu t = 0.75 • • • 3 GeV. The second error, ±0.02% estimates the uncertainty induced by 
the hadronic parameters. 

In view of the smallness of the error, one has to worry about higher order corrections. 
In|l(|, the authors have summed up the leading 0(a n ) corrections for the leading logs in the 
short distance part, which leads to an enhancement of ±0.13%. Furthermore they considere 
the leading QCD short distance corrections of order 0(aa s ), which decrease the short distance 
part by —0.03%. Similarly, 0{ot) n effects should be considered in the long distance part. 

Furthermore, note that we employ mz as an ultra-violet cut-off for the short distance 
corrections, according to the general theorems by Sirlin[[l!| on short-distance electroweak 
corrections to semileptonic processes. But this implies that there is an uncertainty of the 
order of a/(2ir) x 0(1) ~ 0.1% in the definition of the radiative correction, because a change 
of the cut-off scheme would induce a change of the result of this order. Note that the error 
of the 0(a) correction which we have determined is of the same order of magnitude as this 
inherent uncertainty. 

Taking into account the higher order corrections and these uncertainties, we use the fol- 
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lowing value in order to extract fx' 

— (if ^(7)) = (1.3 ± 0.2)% (12) 
1 

Using the 1994 particle data group values for the decay rates and the V us matrix element, 
this yields 

f K = (112.4 ± 0.9) MeV (13) 

From a similar analysis, we obtain 

/tt = (92.1 ± 0.1) MeV (14) 

Two comments are in order: 

• /tt is not defined unambigous at the order 0(a). One could decide to absorb part of 
the radiative correction into /„-. Here we have by definition factored out all radiative 
corrections from f w . Holstein[|[] used a different definition which absorbs terms propor- 
tional to \n(m n /m^) into f n . Use of his definition would yield a slightley lower value of 
f n = (91.9 ± 0.1) MeV. 

• We have estimated the model dependence from the matching uncertainty and the 
hadronic parameter dependence, and so the ±0.1 MeV uncertainty is actually a lower 
limit. 

Taking the full size of the model dependent part as uncertainty, a conservative error estimate 
is 

U = (92.1 ±0.3) MeV (15) 

Note that the error on fx is dominated by the error on |V^ S | rather than by model dependence. 

Our number on f n is to be compared with the particle data group value||, which is based 
on[]l(J. Dividing by y/2 in order to conform with our convention, it reads 

f* pp = (92.4 ± 0.07 ± 0.25) MeV (16) 

where the first error ±0.07 comes from the experimental uncertainty on \ V u d\ and the second 
±0.25 from the matching uncertainty. The numbers are compatible within the error bars. Our 
number has a smaller central value due to the improved matching procedure (our matching 
dependence is in fact an order of magnitude smaller.) 

3.2 The ratios R n and R K 

Let us consider next the ratio Rk of the electronic and muonic modes in kaon decay. From 
Fig. |H we obtain 

SR K = -(3.73 ± 0.02 ± 0.03)% + 0{a 2 ) (17) 

where the error given are the maching uncertainty ±0.02% and the uncertainty from the 
hadronic parameters ±0.03%. These errors give a lower limit of the model dependence of the 
result. 

To get a reliable understanding of the model dependence, we have compared our evaluation 
with the leading model-independent contributions and added up the full sizes of the various 
model dependent parts. This leads to an error estimate of ±0.04%. 
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Figure 2: Radiative correction to the ratio Rk, using different choices for the hadronic pa- 
rameters: Standard choice (I, solid) and variations (II, dashed and III, dotted) 



li cut [GeV] 



Furthermore we have determined the contribution to the loop integral from intermediate 
scales vT^T = 0-5 • • • 3 GeV. This also leads to an error estimate of ±0.04%. 



Finally higher order corrections have to be discussed. In[10| the authors sum up the 
leading logs in m e /m^ to all orders in a, leading to an enhancement of +0.05%. Next-to- 
leading corrections are unknown, and so we allow for an additional ±0.01% uncertainty. This 
leads to our final result 

SR K = -(3.73 ±0.04 ± 0.05 ±0.01 )% = -(3.78 ±0.04)% (18) 
O(a) 0(a n ),n>2 



and 

Similarly, we obtain 



R K = (2.472 ± 0.001) x 10~ 5 (19) 
R w = (1.2354 ± 0.0002) x 10" 4 (20) 



4 — Summary 

Performing an improved matching calculation, we have extracted 

U = (92.1 ± 0.3) MeV 

f K = (112.4 ±0.9) MeV (21) 
where we have, by definition, factored out all 0(a) effects from f n and fx- 
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Furthermore we found 



Rtt — 



r(7r ->• eveiri)) 

r(vr -» ^(7)) 
£(g -> ei/ e ( 7 )) 

r(#-/ii/ M ( 7 )) 



(2.472 ± 0.001) x 10~ 5 



(1.2354 ±0.0002) x 10 



,-4 



(22) 



These are model independent predictions in the sense that the full model dependence is 
included in the error bars given. 

R-x and Rk can be predicted reliably within the standard model, with relative uncertainties 
of a few 10 -4 . Thus they allow for very sensitive low energy precision tests of the standard 
model. 
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